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Abstract 

i The shear viscosity rj has been calculated by using the Green-Kubo relation 

in the framework of a partonic transport approach solved at cascade level. We 
compare the numerical results for r\ obtained from the Green-Kubo correlator 

t—\ with the analytical formulas in both the Relaxation Time Approximation 

(RTA) and the Chapman-Enskog approximation (CE). We investigate and 
emphasize the differences between the isotropic and anisotropic cross sections 
and between the massless and massive particles. We show that in the range of 

Q temperature explored in a Heavy Ion collision and for pQCD-like cross section 

the RTA significantly underestimates the viscosity by about a factor of 2-3, 
while a good agreement is found between the CE approximation and Gree- 

t ^ Kubo relation already at first order of approximation. The agreement with 

^ the CE approximation supplies an analytical formula that allows to develop 

kinetic transport theory at fixed shear viscosity to entropy density ratio, r\j s. 
This open the possibility to explore dissipative non-equilibrium evolution of 
the distribution functions vs T-dependent rj/s and particle momenta in the 
dynamics of the Quark-Gluon Plasma created in ultra-relativistic heavy-ion 
collisions. 



1 Introduction 

The dynamics of the matter created in heavy-ion collisions (HIC) at both the RHIC 
program at BNL and the LHC one at CERN has shown that the azimuthal asym- 
metry in momentum space, namely the elliptic flow i>2, is the largest ever seen in 
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HIC [lj. The comparison of the experimental measured t>2 with hydrodynamical 
calculations has suggested that in these HICs an almost perfect fluid with a very 
small shear viscosity to entropy density ratio rj/s has been created [33 0]. Similar 
conclusions has been obtained also by kinetic transport theory [61 d El E] . At the 
level of comparison performed till now all these approaches have shown an agree- 
ment on the evaluation of the viscosity with Atttj/s ~ 1 — 3. On the other hand the 
differential elliptic flow v 2 (pt) has a very little change from the y/s^N = 39 GeV up 
to the last LHC data at ^s^n = 2.76 TeV. It is an open question if this means an 
equal rj/s of the formed plasma between the RHIC and LHC energies or it is the 
result of different initial conditions and possible different non-equilibrium effects. 

Both Hydrodynamical calculation and transport calculations have been widely 
used to study the dynamics of HIC both showing that the elliptic flow depends 
sensitively on the ratio rj/s. 

Kinetic transport theory starts from a microscopic description of the dynamical 
evolution and hence needs the knowledge of cross section and mean fields, if Vlasov- 
like drifting terms are included [28] . 

Hydrodynamics instead is a macroscopic approach with dynamics ideally con- 
strained by the conservation of energy-momentum tensor and currents. Such over- 
simplification is possible only if dissipation is negligible, while in realistic cases an 
expansion in terms of dissipative terms is needed. This leads to relativistic viscous 
hydrodynamics usually developed according to the Israel-Stewart Theory [3j H]. Here 
viscosity represent an extrinsic parameter input to be calculated by the pertinent 
quantum field theory or related effective lagrangian approximation. 

In kinetic transport theory dissipative dynamics due to finite viscosity is instead 
intrinsically included due to the presence of finite cross sections. However usually 
kinetic theory is applied to the study of HIC starting from the microscopic details of 
the fields and cross sections and it is not discussed directly in terms of viscosity of 
the system. The search for the QGP properties, however, have shown that the shear 
viscosity and in particular the viscosity to entropy density ratio rj/s is a key transport 
coefficient that could be very close to the conjectured lower bound limit, rj/s = 1/4tt. 
This has lead more recently to develop a transport approach at fixed rj/s [HJ El [28] 
that allows to have a direct link to the viscous hydrodynamic language. On the 
other hand kinetic theory is a tool than can allow to investigate the non-equilibrium 
and dissipative effect in a wider range of validity respect to hydrodynamics for both 
rj/s and the momenta of the particles. First attempts in this direction have been 
already developed and applied to the study of the QGP dynamics using the simple 
expression rj/s = 4T/(5a tr ), where Ot r is the transport cross section. However such 
an approach ask for a knowledge of the correct relation between the shear viscosity 
rj and temperature, cross section, mass and density. 

In literature there are several methods for the calculation of the shear viscosity, 
the most employed ones are the Relaxation Time Approximation (RTA) |XTJ and 
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the Chapmann-Enskog (CE) method [12J. The first is based on an ansatz for the 
collision integral in the Boltzmann equation and it does not allow to have control 
of the precision of the approximation. The CE approach is instead a variational 
approach that in principle allows to obtain solutions with an arbitrary accuracy 
which depends on the order of approximation used, [19J, at variance with the RTA 
where it is not possible to have control over the degree of accuracy of the method. 
Nonetheless the RTA approach has been even more often used due to its simplicity 
to evaluate the viscosity for both hadronic and partonic matter [T31 El [TSJ [T6] . 

In this work, we use the Green-Kubo relation [17] to extract the shear viscosity 
of a gluon gas from microscopic transport calculations within a parton cascade 
model with elastic two body collisions. The Green-Kubo method in the contest of a 
relativistic plasma have been discussed also by other groups for the hadronic sector 
pTTj [3D] and also in the QCD sector employing the pQCD matrix elements for 2 — ^ 2 
collisions [10] and more recently also for 2 — y 2 and 2 — y 3 collisions [39J . 

Our objective here is however quite different and it is essentially double folded. 
From one hand we want to compare the two main analytical approximation schemes, 
CE and RTA, with the results obtained evaluating the Green-Kubo correlator solving 
it numerically by the kinetic transport equation in a box with the test particle 
method. From the other hand we will show that for all the case of interest the CE 
already at first order is a pretty good approximation to the exact viscosity. Hence 
providing an analytical relation between rj *-y T, a(9),p, M. Such analytical relation 
supplies a way to construct a kinetic theory at fixed rj/s(T) with a much larger 
accuracy respect to first tentatives [HI El EH] , especially for the case of non-isotropic 
cross section and massive quasi-particles both of interest for a realistic description 
of both quark-gluon plasma and hadronic matter. 

In particular, we will show that the RTA approximation usually employed to eval- 
uate the viscosity in NJL or quasi-particle models [151 12H [25J |2TJ , rj/s = 4T/(5a tr ), 
can lead to large inaccuracy in the r\ evaluation of more than a factor 2, while 
the CE, already at first order, shows a satisfying agreement with the Green-Kubo 
results. 

The paper is organized as follows. In Section II, we discuss the method for cal- 
culating the Green-Kubo correlator paying particular attention to the problem of 
convergency. In Section III, we give an overview of Chapmann-Enskog and relax- 
ation time approximation discussing the comparison with the Green-Kubo method 
for the cases of isotropic and anisotropic cross section and for massive and massless 
particles. In Section IV, we apply our method to the specific case of a gluon plasma 
with pQCD-like cross section that have a Debye screening mass in the gluon propa- 
gator and we show that also in this physical case the relaxation time approximation 
can significantly underestimate the rj/s of the plasma. Finally Section V contains 
summary and conclusions. 
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2 Shear viscosity from Green-Kubo relation 



The transport coefficients like heat-conductivity, bulk and shear viscosity can be 
related to the correlation functions of the corresponding flux or tensor in thermal 
equilibrium |T7]. The underlying physical reason is that dissipation of fluctuations 
have the same physical origin as the relaxation towards equilibrium, hence both dis- 
sipation and relaxation time are determined by the same transport coefficients. Here 
we are interested to the shear viscosity rj for which the Green-Kubo formula assume 
the following expression [18] : 

v = ^J™dtj d 3 x (tt^(x, 0^(0, t)) (1) 

where T is the temperature, n xy is the xy component of the shear component of 
the energy momentum tensor while (...) denotes the ensemble average. In this 
work we determine numerically the correlation function {ir xy (t)7r xy (0)) solving the 
ultra-relativistic Boltzmann transport equation. 

The correlations of shear components can be computed by employing transport 
simulations for a particle system in a static box of volume V at equilibrium. We 
perform such simulation by mean of a relativistic transport code already devel- 
oped to perform studies of the dynamics of heavy-ion collisions at both RHIC and 
LHC energies [HI El EH EH]- The parton cascade developed solves a relativistic 
Boltzmann- Vlasov equation: 



p fi d fl f(x, P )+M(x)d t ,M(x)dgf(x, P )=C(x, P ) (2) 

where f(x,p) is the distribution function for on-shell particles and C(x,p) is the 
Boltzmann-like collision integral that for a one component system can be written in 
a compact way as, 

C(x,p) = J I f{f v f 2l - /i/ 2 )|A^i'2'^i 2 | 2 5 4 (pi +p 2 -p[ -p' 2 ) (3) 

2 1'2' 

where j. = J d 3 P j/(27r) 3 2Ej, M. denotes the transition matrix for the elastic pro- 
cesses and fj are the particle distribution functions, directly linked to the differential 
cross section \Ai\ 2 = lQix s (s — 4M 2 )da/dt with s the Mandelstam invariant. Indeed 
the calculations of the transport coefficients, based on the formal correspondence be- 
tween the kinetic theory and the hydrodynamics, are carried out as approximations 
of the collision integral in Eq.g, see ref.s pH EQl E21 E3 EI] . 

Our aim here is to solve numerically the full collision integral to evaluate the 
viscosity through the Gree-Kubo formula and compare it with the results of the 
RTA and CE approximation scheme. The particle dynamics is simulated via Monte 
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Carlo methods based on the stochastic interpretation of transition JTUJ [8] , according 
to which collision probability of two particle in a cell of volume AV ce u and within a 
time step At is 

p = v "'N^ t Mr; W 

where a to t is the total cross section and v re \ = yjs(s — AM 2 )/2E 1 E 2 denotes the 
relative velocity of the two incoming particle. This method has been shown to be 
appropriate also for collisions between ultra-relativistic massless particles, once one 
choose a sufficiently small time step At respect to the time scale of the process to be 
studied, in our case we have chosen At = 0.01 fm/c. We have performed calculation 
in a stationary box of volume V^ ox = 27 fm 3 with periodic boundary condition and 
a grid cell volume AV ce u = (0.1) 3 fm 3 , checking that further increase of Vb ox and 
reduction of AV ce u does not change the results. The number of test particle N tes t 
allows to reduce spurious statistical fluctuations, usually it is chosen to reproduce 
the behavior of known quantities such as the energy density e and the pressure P. 
We have of course checked that these quantities are correctly reproduced, however 
in our case it is need a much more careful study of the convergency toward the exact 
solution because we are interested to study temporal correlations and not simply 
global thermodynamical observable, see below for a more detailed analysis of this 
point. 

The shear component of the energy momentum tensor is given by 

7T^(x, t) = T^(x, t) = j ^~f{^ P; t) (5) 

where we notice that at equilibrium the shear stress tensor is given by the energy- 
momentum tensor. In our calculation the particles are distributed uniformly in the 
box. Therefore for an homogeneous system the volume averaged shear tensor can 
be written as 

K xy (t)=^it pj f ( 6 ) 

v i=l ^ 

the sum is over all the particles in the box. 

Numerically it is useful to write the correlation function (7r xy (t)7r xy '(0)) as follows 

(7r xy (t)ir xy (0)} = ( lim / dt' Tt xy (t + t')-K xy (t')) = (7) 

\T max ^OG l max Jq j 

Yl ^(iAt + jAt)7r xy (jAt) 



where T max is the maximum time to choose in our simulation, Nx max = T max /At is 
the maximum number of time steps and zAt = t, while (...) denotes the average 
over events generated numerically. 
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Figure 1: Correlation function (n xy (t)n xy (0)) as a function of time for different value 
of the total cross section cr tot , the temperature is fixed to T = 0.4 GeV. These results 
are for massless particles and for isotropic cross section. 

In Fig. (jlj it is shown the correlation function (n xy (t)n xy '(0)) as a function of time 
for different values of the total cross section a tot- In these calculations the particles 
are distributed uniformly in the box in coordinate space, while in the momentum 
space we have considered a thermal distribution (/ ~ e~ E ^ T ) (at a temperature of 
T = 0.4 GeV). In the particular case shown in Fig.Q the particles are massless and 
they interact with an isotropic energy independent differential cross section. As we 
can see the function (ir xy (t)7r xy '(0)) is an exponential decreasing function. We notice 
that the exponential decay with time is not an assumption, but it comes directly 
from the calculation from the dynamical evolution in the simulation. We use this 
fact to fit the correlation function with the following expression, as done in several 



r is the so called relaxation time. As we can see in Fig. ([I]), the initial point of 
the correlation function is independent on the strengh of interaction but it depends 
only on the temperature and volume of the system. This result found numerically is 
indeed expected because it is possible to show that the initial value of the correlation 
function can be expressed by the following relation 



from which we can see that the correlator at t = depends only on thermodynamical 
variables and not on the cross section. We have found an agreement between the 
analytical value in Eq.(|9| and the numerical result at the level of 1%, see Fig. (|2]). 

The slope of the correlation function (i.e. r) depends on the scattering cross 
section and we will show that as expected from kinetic theory r — cj (op) with the 



other works [39j 001 El ED] , 



(n xy (t)iT xy (0)) = (7r xy (0)7r xy (0))e- t/r 





(9) 
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Figure 2: Left: relaxation time r as a function of the maximum time T max . Right: 
initial value of the correlator (TT xy (0)ir xy (0)) as a function of T max . The total cross 
section is fixed to a to t — 0.1 fm 2 and the temperature is T = 0.4 GeV. The dashed 
line is the analytical result. 

coefficient c that depends on the mass of the particles (or better on the m/T ratio, 
see next Section) and on the angular dependence of the scattering cross section. 
Substituting Eq. ^ into Eq. (fll the formula for the shear viscosity becomes: 



This is the formula that we will use in our calculation to extract the shear viscosity. 
The relaxation time r is calculated performing a fit on the temporal range where the 
correlation function assume the exponential form, because at t » t the correlation 
becomes too weak and the fluctuations starts to dominate. A key point is the 
evaluation of the error on the value of the viscosity as coming from the error on the 
initial value of the correlator and the error on the relaxation time r extracted from 
the fit of the correlation function, hence possible deviation from the exponential law 
are evaluated through the error bars themselves. 

In Fig. (|2| it is shown an example of study of the convergence of the relaxation 
time t and the initial value of the correlator (7r xy (0)7r 3:?/ (0)) as a function of the 
maximum time of the simulation T max for the case of isotropic and constant total 
cross section of a to t — 0.1 fm 2 and for a temperature of T = 0.4 GeV. We have 
performed such analysis using a large number of test particles N tes t = 1000, the 
convergency with N tes t will be described soon after. As shown in Fig. [2] (right) 
(n xy (0)n xy '(0)) converges to the analytical value given by Eq.(|9]) (dashed line) with 
good accuracy only at T max « 100 fm/c , which is a very large time scale compared 
respect to the relaxation time for this case, r ~ lfm/c. The r itself is instead 
always quite close to the exact value, but for small T max , however, the evaluated 
error bars show that there is a large uncertainty on the exponential fit that again 
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Figure 3: Left: relaxation time r as a function of the maximum time N tes t- Right: 
initial value of the correlator (7r xy (0)ir xy '(0)) as a function of N test . The total cross 
section is fixed to a to t = 0.1 fm 2 and the temperature is T = 0.4 GeV. The dashed 
line is the analytical result. 

is reduced increasing the maximum time over which the temporal correlations are 
followed. The large error bars are essentially indicating that for small T max it is not 
possible to have a defined exponential decay behavior according to Eq.(|8|. 

We have checked also the convergency vs the number of test particles, in Fig. ^ 
we see that for r it is important to have a large number of test particle N tes t > 500 
to reduce the uncertainty on r at the level of 2%. On the right panel of Fig|3] 
the (ir xy (0)7r xy '(0)) is shown to be nearly indipendent on N tes t, as one can expect 
considering it is the initial value of the correlation hence less affected by the accuracy 
of the system evolution. At variance with the dependence of r on T max the N test 
is of crucial importance. This can be understood because a small number of N tes t 
does not allow to properly map the phase space and consequently the dynamical 
evolution of the system. 

A similar study have been performed for all the numerical calculations shown in 
the paper, we have used T max = 100 r, N test = 1000 and At = 0.01 fm/c. 

3 Chapman-Enskog and Relaxation Time Approx- 
imation vs Green-Kubo 

In literature there are several methods for the calculation of the shear viscosity, 
one of these is the Chapman-Enskog (CE) approach [Til [I~2l [18] and another is the 
relaxation time approximation (RTA) [211 E2]. The difference between these two 
methods resides in the different way in which the collision integral is approximated. 
In the RTA is not possible to have control over the degree of accuracy of the method 
and furthermore it is not possible to go to higher order. In the CE approximation 
instead is possible to obtain solutions with an arbitrary accuracy which depends on 
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the order of approximation used. We show briefly the two different approximation 
methods and formulas for the viscosity that we will use to compare the numerical 
results obtained using the Green-Kubo relation. 

The starting hypothesis of the RTA is that the collision integral can be approx- 
imated by 

c\f] = - J -^— (11) 

r 

where r is the so called relaxation time. In this approximation we are assuming that 
the effect of the collision between the particles is to keep locally the system close to 
the equilibrium. Collisions are assumed to restore the local equilibrium always with 
a relaxation time r that is of the order of the time between two collisions. Hence the 
the distribution function / tends to the equilibrium value f eq exponentially with a 
time t. In this approach it has been demonstrated that the shear viscosity assumes 
the following expression [2H [22] : 



d3 Pa \Pa\ 4 1 req 

15T J (2tt)3 El w a (E a ) Ja 1 ) 

where w a (E a ) is the so called collision frequency and f^ q is the equilibrium distri- 
bution function of particles a with momenta p a and energy E a . The relaxation time 
r a = u(Ea)- 1 is given by 

r~' = M E.) = E \ J ^ §fs g| na, *) C , (13) 

where the quantity W(a, b\c, d) is defined as 

(27r) 4 £ 4 (p a +p b - Pc - 
2E a 2Ef } 2E c 2Ed 



W(a,b\c,d) = r ' v ;;i r :° " \M\ 2 . (14) 



| | 2 is the squared transition amplitude for the 2-body reaction a + b — > c + d. 
The collision frequency Eq.(13) can easily be expressed in terms of the total cross 
section <r tot : 



ME a ) = J j~y 3 2Eam f b a tot , (15) 

that for a constant (energy independent) total cross section becomes energy 
independent and coincides with the standard mean relaxation time 
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T a 1 = W a (E a ) = p a to t (Vrel) (16) 

where (v re i) is the thermal average of the relative velocity and for massless particles 
( v rei) — 1- We notice that in the RTA the interaction appears in the collision fre- 
quency only through the total cross section. However on general physical argument 
the viscosity is expected to depend also on the momentum transfer that on aver- 
age the collisions are able to produce. Therefore one would expect that different 
angular dependent cross sections generate a different viscosity even if they have an 
equal total cross section that is associated only to the probability to collide and 
not to momentum transfer occuring in the collisions. To occur for the q 2 — transfer 
efficiency of the cross section, in transport theory one defines the transport cross 
section, a tT = J d(cos9)a(Q)(l — cos 2 Q), i.e. the differential cross section weighted 
by the momentum transfer that is proportional to (1 — cos 2 0). In the literature 
sometimes to take into account this fact the relaxation time is approximated by 
T tr 1 = (p a tr v rei) j i- e - substituting the total with the transport cross section. This is 



not really coming from the RTA as in Eq.s (12) and (15), however we will also refer 
to it as modified RTA in the following, but we will see that also such an extension 
of the RTA can reasonably approximate the correct viscosity only for the case of 
isotropic cross section, in which case a tr = (2/3)cr tot . 

To describe the CE approach we use the formalism recently developed in Ref. 
|19j . where it is discussed how the shear viscosity at first order [t] s }i can be written 
for the most general case of relativistic particles at finite mass colliding with a 
non-isotropic and energy-dependent cross section as 

T <-y 2 

[Vsh = 777- (17) 
10 coo 

where 70 = — lOh, with z = m/T and h = K^{z) j Ki{z) is the enthalpy and 



rou U> ( W f - \*>? + ) (IS) 



(s) 

The wl are the so-called relativistic omega integrals which are defined as 

/oo 
dyy\y-l) 7 K,{2zy) 



,0) _ / A„,„,U„, 1\7 



K 2 (z) 2 



^2 

X 



f'TT 

I d® sin6<7(s,e)(l -cos* 6) . (19) 
Jo 

where a(s, 0) is the differential cross section and j = | + | (—1)* while y = yfs/IM. 

Note that the angular integration in the u\ for s = 2 is proportional to the transport 
cross section a tr = j dQa(s,Q) sin 2 6. 
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In the following we will consider the case of isotropic cr for both massless and 
massive particles and the more realistic case of non-isotropic cross section for mass- 
less particles. 



3.1 Isotropic cross section: massless case 

In this section we compare our results for 77 using the Green-Kubo relation and the 
analytical result obtained in the relaxation time approximation and the Chapman- 
Enskog approach, for the case of massless gluons r] colliding with an isotropic cross 
section a(s, G) = ctq = const. The CE first order approximation for the shear 



viscosity, Eq.(17), for the case considered is simply: 



IvYce = 0.8— = 1.2— (20) 

O'tr O tot 

where for the last equation we have used a tr = 2/3&tot valid only for isotropic 

cross section. We notice that in CE already at first order of approximation appears 

the differential cross section weighted with (1 — cos 2 G), i.e. the transport cross 

(2) 

section. This comes from the u\ integrals that independently on the i— value are 
proportional to (T tr (s), hence one can re-write coo = f(z)a tr , with f{z) that can 
be written in terms of a combination of Bessel functions with coefficients coming 



from the energy integration present in Eq.(19). So one finds in general n ~ \ jo tr 
as expected. In the next section we derive a more explicit formula for the case of a 
pQCD cross section with HTL-like dressed propagator. 

For this simple case of isotropic cross section and massless particle in the litera- 
ture there exist also higher order calculation up to the most recent work in Ref. [T9] 
where the calculation was extended up to the 16 t/l order. They found that higher 
order approximation converge to the value 

Irftf* = 0.845— = 1.267— (21) 

C tr & tot 

where however approximations beyond the third order give corrections much smaller 
than 1%, while the difference between the I s * and the Vo th order is about 6%. 

In the RTA the shear viscosity for the massless case can be calculated using 



Eq.(12) and Eq.(16) easily obtaining: 



Vrta = 0.8— (22) 

Otot 

Hence the difference between the Chapman-Enskog (CE) at first order and the re- 
laxation time approximation (RTA) can be seen as the substitution of a tot with ot r - 



1 In our result in this Section there is nothing specific of a gluon system and the results discussed 
are however very general. 
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Figure 4: Shear viscosity r\ for a massless system as a function of the total isotropic 
cross section a and for different temperatures. The circles are the results from the 
Green-Kubo method, while the lines are the results obtained using the Chapman- 
Enskog approximation in Eq.(21). 



This result predict a viscosity about 40% smaller respect to the CE approximation. 
As mentioned in the previous section the RTA does not take into account the trans- 
port cross section, but quite often the RTA is used choosing a r _1 = (pcr tr v re i), 
where again for an energy and angular independent cross section one obtains the 
same results in first order CE: 

V*rta = 0.8— = 1.2—. (23) 

"tr O tot 

where we put a "*" to indicate that it is not exactly the RTA but a sort of modified 
one. We will show that this is correct only if the cross section does not depend 
neither on the energy nor on the scattering angle. In other terms it cannot be 
applied to the case of a QCD matter and more generally to hadronic and partonic 
effective lagrangian approach that always predict quite non-isotropic scattering cross 
section. 

Our results obtained using the Green-Kubo formula are shown in Fig. Q by full 



circles and compared to the prediction of CE at 16 th order, Eq.(21 ) shown by different 
lines for each temperature from 0.2 to 0.6 GeV . The error bars for the Green-Kubo 
calculation are small and within the symbols. As we can see we have a very good 
agreement with the analytical results in all the examined range of cross sections 
and temperatures. Discrepancies are observed only at the level of about 2%, which 
is comparable to the uncertainty in the numerical evaluation of the Green-Kubo 
correlator. 

In our simulations the gluons are considered as massless particles and due to the 
fact that they are distributed according to the Boltzmann distribution the entropy 
density is given by s = 4p = 4d g T 3 /ir 2 with d g = 16. We have checked that the 
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Figure 5: rj/s as a function of the temperature and for different cross section a tot = 
1, 3, 6, 9, 30 mb. The circles are the results of the Green-Kubo method while the lines 
are the result obtained using Eq.(24) from CE approximation. 



transport code correctly reproduce the analytical formula by calculating (II 00 + 
U xx )/T. Therefore for a temperature independent cross section we have 

rj/s = 0.195^— (24) 

In Fig. ([5]) the shear viscosity to entropy density ratio rj/s is shown as a function 
of the temperature and for different total cross section a. The circles are the results 



of our simulations while the lines are the analytical results obtained using Eq.(24). 
From the value of rj/s we can see that we have explored the scale of a, p, T of 
interest for the study of the QGP. The viscosity extracted with the Green-Kubo 
relation is in good agreement with the CE approximation in a very wide range of 
density and cross section. We notice that the results shown correspond to a gas 
condition with a mean free path 0.005 < A < 5.95 fm corresponding respectively 
to a — 30 mb and p = 45.6 fm -3 (T = 0.6 GeV) and a = lmb and p = 1.68 fm -3 
(T = 0.2 GeV). Sometimes in the literature it is stated that the Boltzmann transport 
equation is applicable only in a regime of diluteness. However the meaning of such 
a statements is often misunderstood. It does not mean that one cannot study by 
mean of Boltzmann transport equation the system dynamics if the mean free path 
is small respect to the average particle distance. In fact our work shows how one 
can fix the cross section according to the CE approximation and the Green-Kubo 
formula, demonstrating that in such a way one really has a fluid with the wanted 
viscosity even if this corresponds to very small mean free path respect to particle 
distance. Of course what one cannot deduce from this is if microscopically in the 
real system the viscosity comes only from two-body collisions. 
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Figure 6: Left: shear viscosity t] as a function of the total cross section for different 
value of the mass m = 0,T, 5T, 10T and for a fixed temperature T = 0.4GeV. 
The circles are the results obtained using the Green-Kubo relation while the lines 
are the analytical results obtained using the Chapman-Enskog approximation to 
the first order of approximation Eq.(17). Right: f(z) = T)a to t/T as a function of 
z = m/T, the diamonds are the results obtained using the Green-Kubo relation 
while the solid and dashed lines are given by the CE approximation respectively for 



the first (dashed line) and second order (solid line). The modified RTA Eq.s(12) 
with transport relaxation time is shown by dashed dotted line and standard RTA 
by the dotted line. 



3.2 Isotropic cross section: massive case 

The calculations discussed till now consider only massless particles: an interesting 
extension is that with finite mass. The importance of this kind of calculation for a 
partonic system is that it can be employed to evaluate the viscosity in quasi-particle 
models [331 EI] • We remind that the quasi-particle approach is important because 
it allows to describe the equation of state obtained on the lattice of a plasma of 
quarks and gluons [281 EE]. This is achieved incorporating the non-perturbative 
effects of QCD in a temperature dependent mass term for both gluons and quarks. 
Furthermore finite mass is relevant also if one wants to study the viscosity of an 
hadronic fluid or employ transport theory to study the dynamics of hadronic matter 
the mass cannot be neglected. 

In the left panel of Fig.(|6]) the results for the shear viscosity rj for the massive case, 
obtained using the Green-Kubo relation, are shown and compared to the analytical 
result obtained using the CE approximation. For a constant differential cross section 
a(s, 0) = o"o the shear viscosity in the CE approach can be written at first order 



approximation, Eq.(17), in terms of the total cross section a tot as 



T 

<Jtot 



(25) 
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where the function f(z) [T2] 120] g| 



/(*) 
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16 (15z 2 + 2)ir 2 (2^) + (3z 3 + A9z)K 3 (2z) 



(26) 



where f(z) — >■ 1.2 in the ultra-relativistic limit z — >■ 0. In the right panel /(z) of 
Fig. ^ is shown as a function of z = m/T for the first order (dashed line) and second 
order (solid line) of approximation. The diamonds represent the mean value over 
the different cross sections of the results obtained using the Green-Kubo relation for 
the viscosity 77 for each value of mass explored. As we can see, we have a fairly good 
agreement with the analytical result already at first order. However as already seen 
for M = there is a discrepancy of about 6% at first order corresponding to the 
difference between r)a to t/T = 1.2 and 1.267. Here we have evaluated numerically 
also the CE expansion at second order that according to [19J: 



T 7q en - 2 70 71 coi + 7? c 00 



10 



C 00 C 11 ~~ c i 



01 



(27) 



where 
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44 
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175 



32 



2^0 



(2) 
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The behavior of [^ s ] 7/ cr tot /T is shown in Fig. ^ (right) by solid line. We see 
that the second order CE improves the agreement with the Green-Kubo results 



2 In the De Groot et al. book it is present a typos in the denominator where 15 z 2 is written as 



5z A 
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reducing the discrepancies to less than 3%. It seems that with increasing mass the 
convergency is slower and higher order are more relevant. 

Our result shows that for massive particle the coefficient between t] and T/a to t 
increases by about 40% for value of M/T typical of gluon mass in quasi-particle 
models [331 H22 ESI EZ] or protons around and below the cross-over temperature. 
The comparison with the RTA shows also in this case that as we move from the 
simplest case of massless particle the viscosity is underestimated. In the case of 
the pure RTA (dotted line) the correct value, about a 40% larger, become about a 
70% larger for a mass M = 10T, but also the modified RTA that is equal to CE 
at first order for M = under predict the r] when a finite mass is considered. The 
discrepancy increase with the mass and it is about 15% for M = 10T even if the 
cross section is isotropic. 



3.3 Anisotropic cross section 

In this section our aim is to study the more realistic case of angular dependent cross 
section, being in the context of a gluon plasma we choose typical elastic pQCD 
inspired cross section with the infrared singularity regularized by Debye thermal 
mass rriE>: 

da _ 97rag 1 / m 2 D \ 

dt~ 2 (t-mtf{ 1+ s ) (31) 

where s, t are the Mandelstam variables. Such kind of cross sections are those 
typically used in parton cascade approaches [321 H21 El |9l [271 1101 I7J- The total 



cross section corresponding to Eq. (31) is a to t = 9vro;^/(2m| ) ) which is energy and 
temperature independent. Here our objective is not to estimate the value of rj (see 
next Section) but only to explore the impact of non-isotropic cross section on the 



comparison among RTA, CE and Green-Kubo methods. In Eq.(31) the Debye mass 
rriu is a parameter that regulates the anisotropy of the scattering cross section. We 
will vary it to regulate the anisotropy, but fixing the total cross section constant by 
keeping constant the a 2 s /m 2 D ratio. We note that for a plasma at temperature T 
the average possible momentum transfer q 2 w (3T) 2 , hence for mo » 3T Eq. (31) 



acts as an almost isotropic cross section and we should recover the results of the 
previous Section. On the other hand, we notice that the well known HTL estimate 
of a gluon plasma viscosity [H] is valid only in the limit of g = mn/T « 1, i.e. for 
very anisotropic cross section. 

We have seen that the modified RTA gives the same result of the first order CE 
for isotropic cross section. In order to perform the same analysis for the non-isotropic 
case we have to calculate the transport cross section a tr : 

/dcj 
— sin 2 6 dt = a to t h(a) (32) 
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where h(a) = 4a(l + a) [(2a + l)Zn(l + 1/a) — 2] and a = m 2 D /s. For — >■ oo 
the function h(a) — >■ 2/3 and we recover the isotropic limit, <r tr . = (2/3)cr toi , while 
for finite value of mjj the function h(a) < 2/3. One can generalize the concept of 
transport relaxation time r^ 1 = (pcrtrVrei) that for isotropic cross section is equal 
to (2/3)(pa to t Vrei) writing 

T tr l = (p °tr V r ei) = P CTfot (v r el h(d)) (33) 

The shear viscosity t]* RTA in the modified RTA is therefore given by 

Vrta = 0-8 7 ]-- — (34) 

{v re ih(a)} a tot 

This is the same formula used in several transport calculation to fix the viscosity 
(321 El EH], but also the formula used to evaluate the viscosity in effective lagrangian 
approaches [13] . 

In the CE approximation the situation is more complicate and it is not possible 
to write the viscosity in a similar way because the differential cross section enter in 



the relativistic omega integrals oj\ s \ Eq.(19), where the function h(a) is multiplied 



by integrated Bessel functions. However one can rewrite Eq.(17) in the massless 
limit (z — > 0) in the following form 

l Vs y cE = 0.8g(a,T)— (35) 

Ptot 

of course if g(a, T) = {v Td /i(a)) -1 the CE at first order and the RTA would give the 
same prediction. However, as can be seen from the expression of coo, the g(a,T) is 
a much more complex expression of which the function h(a), relating the transport 
to the total cross section, is only a small part. 

In the left panel of Fig. ^ it is shown the shear viscosity r] as a function of 
the Debye mass at fixed total cross section at a to t = 3mb and for three different 
temperatures T = 0.3 GeV, T = 0.4 GeV and T = 0.5 GeV. The solid, dashed and 



dot-dashed lines are the behavior of rj in the modified RTA approximation, Eq.(34), 
while the symbols are the result with the Green-Kubo formula. It is evident that 
there is a strong disagreement between the two as soon as nip is such to move from 
the isotropic limit indicated to guide the eye by dotted lines (which corresponds to 
the standard RTA). Therefore, we see in general that even if the total cross section 
is kept constant, the anisotropy of the cross section cause a strong enhancement 
of the viscosity rj. However the increase is very strong and the difference between 
transport and total cross section is not able to account fully for such an increase. 
On the right panel of Fig. ^ we compare the Green-Kubo results (symbols) with 
the prediction of CE at first order (solid, dashed and dot-dashed lines). In this case 
we find a very good agreement between the two, hence the CE already at first order 
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Figure 7: Left: Shear viscosity r] as a function of the Debye mass mo for two 
different values of the temperature T = 0.3 GeV (blue thick lines) and T = OAGeV 
(red thin lines). The dotted line is the isotropic limit when ^> 1 GeV, the solid 
line is the relaxation time approximation with r _1 = (patrVrei) and dashed line 
is the Chap man- Enskog approximation at firt order. Right: Dashed line represent 
Chapman-Enskog approximation at the first order for three different temperature 
(T = 0.3, 0.4, 0.5 GeV) while the open circles are the results obtained using the 
Green-Kubo relation. 



is able to account for the correct value of r] even if the cross section is so forward- 
backward peaked to cause an increase of about an order of magnitude respect to 
the same total cross section but isotropic. The RTA approximation would severely 
underestimate the viscosity. We can also see that for mjj ~ 8 — 10T the isotropic 
limit is recovered and both CE and RTA coincide but this is essentially the limit 
discussed in the previous Section. We note that the calculation have been performed 
down to quite low value of the screening mass, = 0.1 GeV. This for T = 0.5 GeV 
would correspond to anisotropic cross section that in the HTL approach corresponds 
to g = ttid/T = 0.2. Nonetheless within a precision of about 5% the first order CE 
is able to account for the correct value of r] even for such forward-bacward peaked 
cross section. This result further validates the approach in Ref . |44j . 



4 Viscosity of a Gluon Plasma 

In this last part of the paper we want to apply the analysis performed for quite gen- 
eral case to a realistic case of the shear viscosity to entropy density ratio for a gluon 
plasma. Therefore we consider massless gluons in thermal equilibrium interacting 
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via two-body collisions corresponding to direct u— and t— channel: 



da 99 ^ 99 _ 97ra 2 1 



dq 2 2 (q 2 + m 2 D 



where mp is the Debye mass, rrin = Ty/Aira s according to HTL calculations and 



«s(T)= — 2 , A = 200 Me V (37) 



is the pQCD running coupling costant. A full pQCD calculation with HTL dressed 
propagator would include also the s— channel and all the interference terms, however 
it has been shown that the u— and t— channels are the dominant ones [44]. On the 
other a full HTL calculation as in Ref. [44] would require an energy dependent propa- 
gator with both longitudinal and transverse components, however our objective here 
it is not to have the best evaluation of rj/s, but to discuss the comparison between 
CE, RTA and the numerical Green-Kubo method for a quite realistic case and also 
under the same condition of previous work based on the Green-Kubo method [40] . 
The total cross section in this scheme is energy and temperature dependent: 

_ 9na 2 s s 
2 nip s + m D 

For this realistic case we have calculated the rj/s as a function of the temperature 
for different temperatures in the range 0.2 <T< 1.0 GeV by mean of the numerical 
Green-Kubo method. The results shown by red symbols in Fig. [8] are compared 
to the CE (red dashed line), the modified RTA (red dashed dot line) described 
above and the simple RTA ( black double dash-dot line) with the relaxation time 
proportional to the total cross section and not to the transport one. We find again 
that the CE again is in excellent agreement with the Green-Kubo result at the 
level of accuracy of 4%, while the modified RTA significantly underestimates the 
rj/s by about a 20% at T ~ 0.2 GeV. At increasing temperature the discrepancy 
tends to increase up to about a 50%; we can understand this result on the base of 
the discussion in the previous section because the m 2 D /T 2 cx a s (T) ratio becomes 
smaller at increasing temperature and therefore the cross section appear effectively 
less anisotropic. 

We notice that for T ~ 0.2 GeV we predict a rj/s ~ 1 in approximate agreement 
with the extrapolation of calculations in Ref. [44], even if however we have not used 
the full HTL propagator. 

A significant discrepancy is found respect to Ref. |40] especially for small tem- 
peratures, while for temperatures above T > 0.5 GeV there is a tendency toward an 
agreement. We have however used the same cross section in [40] hence in principle 
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Figure 8: Shear viscosity to entropy density ratio rj/s of a gluon plasma interacting 
through the differential cross section in Eq. (36) as a function of temperature T: 



red circles are the results obtained using Green-Kubo correlator, red dashed line 
represents Chapman-Enskog first order approximation. The dot-dashed line is the 
prediction with the modified RTA approximation and the dotted lines is the standard 
RTA. 



we should have a very similar r]/s(T). However for T < 0.5 GeV they found an in- 
creasing of r]/s, while our results have a monotonic logarithm increase as a function 
of temperature as one would expect from the running coupling a s (T) and the T 
dependence of the screening mass. It appears difficult to justify even qualitatively 
a non monotonic behavior of r)/s(T). 



5 Conclusions 

The recent developments in the physics of strong interaction at finite temperature 
T and baryon chemical potential fis is focusing many efforts on the calculation of 
transport coefficients and in particular of the shear viscosity There are several ap- 
proximation schemes that relates the viscosity rj to the microscopic cross section. 
The validity of such approximations, at least at the typical scales of quark-gluon 
plasma, was not yet investigated. We have developed a method to solve numerically 
the Green-Kubo formula for the shear viscosity for the case of a relativistic Boltz- 
mann gas. Our objective was to compare the Chapman-Enskog approximation and 
the relaxation time approximation with the result from the Green-Kubo correlator 
that in principle should give the correct result. We have performed such a compari- 
son from a variety of physical case: isotropic and non-isotropic cross section massless 
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and massive particles at different temperatures. Our work shows that the relaxation 
time approximation with a relaxation time r oc a tot always underestimate the shear 
viscosity even by more than a factor of 2-3. Moreover even the modified RTA with 
t oc o tr gives a satisfying prediction only in the unrealistic case of massless particles 
and isotropic cross section. For realistic pQCD-like cross section one has on average 
a 30 — 40% underestimate for rj respect to the Green-Kubo calculations. Instead the 
Chapman-Enskog appears to be a much better approximation scheme that already 
at first order have shown an agreement at the level of 4% with the numerical calcu- 
lations of the Green-Kubo correlator for all the physical case considered including 
very forward peaked cross section and massive particles. 

We note that the agreement of CE approximation at first order with the Green- 
Kubo method also supplies a relatively simple analytical expression that can allow 
to developed kinetic transport theory at fixed viscosity with very good precision. 
Our work shows that the current used simple relation rj/s = 0.8T/<7j r can lead 
to significantly underestimate, even by about a factor of two, the rj/s. Therefore 
previous works (HI EHl [26] in this direction are only approximately valid. 

In the next future it would be interesting to extend the study also to a collision 
integral including three-body collisions whose role in the equilibration dynamics in 
ultra-relativistic heavy-ion collisions [TUJ, [TJ H5] and in the determination of the rj/s 
is of current interest. 

Finally, it would be interesting to apply the method for the calculation of rj/s 
for a quark-gluon plasma, this could allow also a more precise evaluation of this 
quantity in quasi-particle models where till now only rough approximation for the 
relaxation time have been considered [H)J 1221 EHl US] • 
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